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It is shown that the criticism presented in the Comment by Galanakis et al on the paper by 
Efetov et al Q is irrelevant to the bosonization approach. 



In spite of the far going conclusion that our bosoniza- 
tion scheme cannot help to overcome the negative sign 
problem, the Comment does not really address the suit- 
ability of the method for Monte Carlo (MC) calculations. 
We wrote in the paper that our mapping of the fermionic 
model onto a bosonic one was exact having in mind 
the continuous with respect to (imaginary) time limit. 
Any quantum MC scheme implies generically a discrete 
time with extremely strong variation of the Hubbard- 
Stratonovich (HS) field 4> on different slices. In this case, 
our bosonization scheme is not necessarily exact for a 
given HS field and we make an approximation. If ev- 
erything remained exact also in this case, we would not 
be able to overcome the sign problem because the parti- 
tion function Z [(p] for a given (j) would be the same for 
both bosonic and fermionic representations and could be 
negative. 

In our opinion, the comment of Galanakis et al. is 
based on a misunderstanding of the definition of Zt- In- 
deed Zh should not be understood as coming from Eqn (9) 
but from Eqns(12-13). [We apologize for this confusion, 
which is the result of writing the paper several times in 
order to improve the presentation]. Having failed to un- 
derstand this point, Galanakis et al. re-derive the same 
steps leading to Zj and jump on the far reaching con- 
clusion that our technique does not solve the MC sign 
problem. However, the bosonization method starts later 
with Eqs. (12, 13), which means that Galanakis et al 
make their conclusions not about the bosonized model 
but still about the original fermionic one. 

They fail as well to understand that for finite time 
slices, our method is not exact but requires an approxi- 
mation. The crucial approximation is made when writing 
Eq. (13). According to Eq. (11), the function Ary (t) 
is expressed in terms of the Green functions at slightly 
different times r and t -|- (5, where S — > +0. Therefore, 
deriving the equation for the function Ary and making 
no approximations one would have the fields (f>r and (pr' 
at slightly different times r and t + 5. Nevertheless, we 
put i5 = in Eq. (13). In the continuous limit, (imply- 
ing subsequent averaging over (pr (t)), this approximation 
becomes exact. Therefore our field theory based on the 
introduction of superfields is exact. 

However, working with finite slices of time and putting 
S — changes the function Zb [(j)] for a given configuration 
of (p. Taking Eq. (13) with the fields (t) and (t) at 



coinciding times r results in a symmetry of the solution 
Ar.r' (t) under the replacement r ^ r' . We argue in the 
paragraph after Eq. (23) that any possible singularity in 
the integral over u in the exponent should be absent due 
to this symmetry. This means that the imaginary part in 
the exponent does not arise and the function Zf, [0] must 
be positive. 

In order to make everything well defined one needs 
a regularization, otherwise the solution of Eq. (13) is 
not uniquely defined. [The fact that the solution of Eq. 
(13) is not uniquely defined is actually remarked in the 
Comment in the paragraph before the conclusion] . In our 
PRL paper we had no possibility to discuss this question 
due to the lack of the space. However, we emphasize that 
Eq. (13) itself and the way we choose its solutions lead 
to a real exponent in Eq. (12) and to a positive Zb [ip], 
thus providing a method of calculations free of the sign 
problem. 

The regularization of Eq. (13) is absolutely necessary 
before this equation can be used for explicit numerical 
computation. This goal is achieved in our recent preprint 
3] where all details of the derivation of Eq. (13) are dis- 
cussed and an explicit regularization is introduced. We 
make a check of our scheme by explicit diagrammatic ex- 
pansions and demonstrate perfect agreement. We give 
detailed explanations about the role of the symmetry 
r ^ r' for removing the singularity of regularized solution 
of Eqn (13). We present there an explicit formula, Eq. 
(4.49), that can directly be used for numerics. We have 
checked it analytically and numerically for a static (p^. 
This formula gives by construction real positive Z^ [(p] for 
any time dependent fields (pr (t) , property which should 
make the computation free of the sign problem. 

We conclude that the claim of the Comment that our 
method cannot resolve the sign problem is not correct. 
Moreover, it is not based on a discussion of the bosoniza- 
tion scheme. The authors of the comment have repeated 
the derivation of the equations preceding the key equa- 
tion (13) of [H where the true bosonization starts. Al- 
though their discussion helps to correct several inaccura- 
cies in the PRL paper, they miss the point that a cru- 
cial approximation is made for the piece- wise fields (pr (t) 
used in the MC calculations and thus come to an incor- 
rect conclusion. As it stands, the Comment of Galanakis 
et al. is not scientifically valid. 
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